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Abstract 

We study the wealth distribution of the Bouchaud-Mezard (BM) model on complex networks. It 
has been known that this distribution depends on the topology of network by numerical simulations, 
however, no one have succeeded to explain it. Using "adiabatic" and "independent" assumptions 
along with the central-limit theorem, we derive equations that determine the probability distribu- 
tion function. The results are compared to those of simulations for various networks. We find good 
agreement between our theory and the simulations, except the case of Watts-Strogatz networks 
with a low rewiring rate, due to the breakdown of independent assumption. 
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I. INTRODUCTION 



It is well known that the toplogy of the network changes the dynamics dramatically. 
After the discovery of the absence of epidemic threshold in a scale-free network[l|, many 
researchers have focused on the dynamics on complex network, such as synchronization 0, [J , 
pattern formation^], and other phenomena. 

In this paper, we study the Bouchaud-Mezard (BM) model on complex network^]. The 
power-law behavior of wealth distribution, Pareto's distribution, has for over a century been 
one of the main cornerstones of econophysics[6[]. One of the simplest models that explain 
this law is that proposed by Bouchaud and Mezard, which consists of multiplicative noise 
and globally coupled diffusion. After the proposal of the BM model, several researchers nu- 
merically investigated the generalized BM model in which diffusion occurs between adjacent 
nodes in a complex network^, Isj]. While, the research revealed that the network topology al- 
ters the wealth distribution, there has been no quantitative theory that sufficiently explains 
these simulation results. 

Recently, we proposed a new theory for the BM model on a random network jjj. Using 
several assumptions that we describe in the next section, we derived equations that determine 
the static probability density function (PDF) of the wealth. The results of this analysis were 
compared with those of the numerical simulations and good agreement was obtained. 

However, the wealth distributions were analyzed only on a random network in the previous 
paper, and the distributions on other complex networks was left as an open problem. The 
aim of this paper is develop our previous work so that it is applicable to a general complex 
network. Using the same techniques applied in the previous paper, we derive the equations 
that determine the static PDF for the BM model on a complex network for a given adjacency 
matrix. The results are evaluated by a comparison with the numerical simulation results, 
and our method is revealed to perform well for many network systems. 

This paper is organized as follows. In section [Til we derive the self-consistent equations 
that determine the PDF of the BM model on a complex network whose adjacency matrix is 
given. The results obtained with this PDF are tested in sec. IHIl by comparing with simula- 
tions on several networks such as a random network, the Watts-Strogatz (WS) network, and 
a real social network. Finally, we discuss the results and future problems and summarize 
the paper in sec. IV. 
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II. THEORY 



We consider the BM model on a undirected complex network consisting of iV nodes, 
whose adjacency matrix is given by A = (ay). The dynamics of x±, i.e., the wealth on node 
i, are determined by the following Ito-type stochastic differential equations. 

N 

dxi = a ij( x j ~~ x i)dt + V2adWi, (1) 
i=i 

where J, a and Wi represents the diffusion constant of wealth, strength of noise, and the 
standard Brownian motion, respectively. This equation has no static PDF; however, the 
normalized wealth xj (x), where (■ ■ •) represents an average over all nodes, can have a static 
PDF. For simplicity, x is used in place of xj (x) in the following discussion. 

To obtain a static Xi distribution, we make the "adiabatic and independent" approxima- 
tion introduced in our previous paper {9]. We first assume that all Xi values are independent 
and that the correlation between different nodes is negligible. Under this assumption, we 
assume the PDF can be decomposed as p(xi, x 2 , • • • , x^) = pi(xi)p 2 (x 2 ) • • • Pn(%n)- Second, 
we assume that the rate of change of Xi = j- J2j dijXj, the average x around node i, is much 
slower than that of Xi, where dj represents the degree of node i. Using these assumptions, 
we can calculate the static PDF as follows. 

Suppose Xi is constant. Then pi(x\xi), the conditional PDF of x on node i for a given Xi, 
is a solution of the following Fokker-Planck equation. 
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at ox ox z 
The static solution of this equation, p eq ,i(x\xi), is given by 



p eq ,i(x\xi) = C di (xi)exp(-a di Xi/x)x , (3) 

where a d = Jd/a 2 and C d (x) = {a d x) l+ad /Y(l + a d ). Using the adiabatic assumption, the 
static PDF of x at node i is given by 



Peq,i{x) = J dxPi(x)p eq4 (x\x) , (4) 

where Pi(x) represents the PDF of X{. 

To obtain P{(x), we use the independent assumption again together with the central- 
limit theorem. Under this assumption, Pi(x) can be assumed to be the distribution of the 



average of the independent variables Xj, where j represents a node adjacent to node i. If the 
variances of Xj are finite for all j, we can apply the central-limit theorem to approximate 
Pi(x) by the Gaussian distribution 



Pi (x) = . — — exp 



{x - M) 2 
2S? 



(5) 



where M, and Si represent the average and variance of Xj, respectively. From Lindeberg's 
theorem, these values are given by 

Mi = ^2a ij {i j /d i , (6) 

j 

and 

S? = E^iK, (7) 
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where \ij and s 2 are the average and variance of Xj, respectively. 

We can determine Mj and Si by calculating /ij and s 2 for x using Eqs. (EJ), Q, ©, and 
(J7J). By using J °° dxp eq> i(x\x)x = x and J °° dxp eq ^{x\x)x 2 = ^ di _ x x 2 for a > 1, we find 

/*oo />oo /*oo /*oo 

/ dxxpi(x) = / dxPi(x) / = / dxxPi(x) ~ Mj , (8) 

JO JO JO JO 

and 

/ dxx 2 Pi {x) = r dx^^-x'Pix) ~ (m + S?) , (9) 

J Jo a di - 1 a di - 1 v / 

where we assume Mj 3> We thus obtain Mj = /ij = 1 for all i, and S 2 is the solution of 

s 2 = 2*5 + 1 . (10) 
®di - 1 

Using Eq. (jTJ), we can rewrite this equation as 



For convenience and in anticipation for later analysis, we define a new matrix R = (r^) as 

Tij = -1)°* ■ (12) 

We then finally obtain the following equation for s 2 : 

£(<% - r^sj = 1/K - 1) . (13) 



Eqs. ©, gD, ©, and (USD define the PDF of x. 

Before concluding this section, we note in the following about wealth condensation, de- 
fined as the divergence of sf. Because s 2 > for all i, Eq. (fT31) must have a solution s| > 
to have any meaningful result. If J is large enough, this condition is satisfied provided that 
min((ij) > 2. For a J value that satisfies > 1 for all i, Rij is a non-negative matrix. 
Therefore, the solution of Eq. (fl3l) satisfies sf > if p(-R) < 1, where p(R) represents the 
spectral radius of R. From Frobenius's theorem, p(R) is equal to the largest eigenvalue of 
R, and thus we conclude that Eq. ([13]) has a meaningful solution if the largest eigenvalue 
of R is less than 1. Noting that di = J2j a ij, p(R) < rnax(J]j an d ry — >■ a^/d 2 for 

J — )■ oo, p(-R) is less than 1 for large J, and a non-condensed phase appears, provided that 
min(<ij) > 2. 

In contrast, we always have a condensed phase at small J due to the divergence of the 
Tij elements. As we decrease J from a large value, the values increase and diverge when 
aft = 1; p(R) also diverges in this case. Therefore, the wealth is condensed if J /a 2 < 
max(l/dj). 



III. SIMULATIONS ON VARIOUS NETWORK MODELS 

To test our theory, we compare the results with numerical simulation results on several 
networks. We begin with the network created by the Erdos-Renyi algorithm and then 
consider several networks created by the WS model[lo| to test the effect of clustering. It 
is shown that our theory fails when the rewiring rate p is too small, but good agreement is 
obtained when p > 0.1. Finally, we test our theory on a real network. It would be ideal to 



apply the theory to a real economic network, but as there is currently no avai 



able data, we 



use the American college football network obtained by Girvan and Newman instead. For 
this case, the PDF obtained from the simulations shows good agreement with our theory. 
For all the simulations, a 2 is set to 1. 



A. Random network 

The network created by the Erdos-Renyi algorithm, shown in Fig. (TJ has 200 nodes and 
a mean degree of (k) = 10. 
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FIG. 1: The random network used for the simulations in sec. III. A 
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FIG. 2: A log- log plot of the PDF on the random network obtained by numerical simulation for 
J = 0.3 (left), 0.5 (middle), and 2.0(right). 

In Fig. |2] we plot the PDFs obtained by numerical simulation for J = 0.3, 0.5, and 
2.0. The PDFs obtained using our theory are indicated by the solid line and show good 
agreement with the numerical simulation for all x for the J = 2.0 case, but there are minor 
discrepancies at small x values in the cases of J = 0.3 and 0.5. The discrepancy is due to 
the Mi ^> assumption used for deriving Eqs. (EJ) and ([9]). The maximum Sf value 
calculated by Eqs. (03) and is 0.33 for J = 0.3 and 0.09 for J = 0.5, while M { = 1. 
Thus, the M, 3> Si approximation is not valid as y/S~i is larger than 0.3. However, the 
approximation is valid for the J = 2.0 the maximum Sf value is 0.01. 




FIG. 3: Upper: The small-world networks created by the Watts-Strogatz algorithm with rewiring 
probabilities of p = 0.01,0.1, and 0.2 used to test our theory. Lower: The PDFs obtained from 
simulations on each networks with J = 2.0. The solid lines show the results of the method proposed 
here. 

B. Small- World Network 

Next, we consider three networks created by the WS algorithm. The networks have 
iV = 100 nodes and a mean degree of (d) = 8 with different rewiring rates of p = 0.01, 0.1, 
and 0.2. Fig. [3] shows both the networks and PDFs. In the case of p = 0.01, our theory 
differs from the numerical simulation, while we find good agreement for the p = 0.1 and 
p = 0.2 cases. 

To investigate the reason for this discrepancy, we examine the scatter plots of x and 
x for p = 0.01 and p = 0.1 in Fig. HJ For p = 0.01, we find that there is a strong 
correlation between x and x that weakens for p = 0.1. The differences can be quantified by 
the correlation, which is 0.74 for p = 0.01 and 0.48 for p = 0.1. This result suggests that 
the discrepancy is caused by the breakdown of independent assumption. 

It is reasonable to question here when and why the independent assumption fails. The 
small-world property does not appear to affect the independence of wealth because the 
networks show high-clustering and small betweenness for both the p = 0.01 and 0.1 cases. 
We return to this point in the discussion section. 
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FIG. 4: Scatter plots of x and x for the BM model on WS networks with p = 0.01 and 0.1 and 
with J = 2. 

C. American college football network 



The American college football network obtained by Girvan and Newman[ll| includes 115 
nodes and 613 edges. The simulated PDFs for J = 0.3 and 0.5 are shown in Fig. |5j In both 
cases, our theory reproduces the PDF well. 

The agreement between theory and simulation seems excellent. We have not understand 
the reason for this good agreement, however, one possible reason is that this network has 
community structure. As Girvan and Newman showed, this network has clear communities. 
As we will see in the later section, the spatial correlation of wealth would be small for 
globally coupled network. Therefore it is natural to assume that wealth correlation is also 
small, if the network is consisted from several densely connected subgraphs. However, we 
still need further study to explain this excellent agreement. 



IV. DISCUSSION AND CONCLUSIONS 



Here, we have developed a theory for the BM model on complex networks. By generalizing 
our previous work, we were able to propose a theory applicable to a complex network with 
a given adjacency matrix. Using the adiabatic and independent assumptions, we derived 
the equations that determine the static PDF of the wealth. The result was compared to 
numerical simulations, and we found our theory works well on a Erdos-Renyi network, WS 
networks with p > 0.1, and the American college football network. We also found that the 
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FIG. 5: The simulated PDFs of the BM model on the American college football network for J = 0.3 
and 0.5. The solid lines represent the results obtained by our theory. 

theory did not perform well for a WS network with p = 0.01 as the independent assumption 
is no longer valid. 

The theory does not apply to a WS network with a small p value because of the large 
spatial correlation in this case. Thus, we may ask when and why the spatial correlation 
becomes large. While it is difficult to give a complete answer to this question, the following 
discussion suggests that the properties of the Laplacian are essential to this problem. 

Suppose that J is very large and the fluctuations of x are very small. Under such cir- 
cumstances, we can approximate xdW in Eq. ([[]) as dW and obtain 

dxi = - J LijXjdt + \f2adW , (14) 

3 

where L = (L^) is the Laplacian matrix L%j = Sijdi — a^. To investigate the correlation 
in this toy model, we assume for simplicity that the network is undirected and connected. 
Introducing the eigenvalues and corresponding normalized eigenvectors of L as Ai = < 
A 2 < •••Ajv and v 1; v 2 , ■ ■ • v^r, we can write (xi,x 2 , ■ ■ • -,x^) T = Y.iLiQi^i- Then L is 
diagonalized and Eq. (fl4l can be written as 

dq { = -J\ t q t dt + V2adW . (15) 

The static distribution of g«(2 < i < N) is calculated as 

(16) 

The qi term does not have a static distribution because Ai = 0. However, because v x = 
-^=(1, 1, • • • , 1), this term gives the uniform change, Xi — > Xi + C, and does not contribute 
to the correlation. 



P{<b 



2,71(7 



exp 



2a 2 
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From Eq. (fi~6|) and the relation X{ = J2j Qj{ w j)ii we obtain the covariance 

N 2 

((*, - ( x ))( Xj - (x))) = £ ^-(yi)i(yi)i • (i7) 

1=2 JX l 

Our independent assumption implies that the off-diagonal elements of Eq. ( !T7|) are smaller 
than the diagonal elements. This condition is satisfied if every Vj is "localized", in other 
words, has only one large components. It has been shown by some researchers that eigen- 
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13l | . Our theory will be applicable to 



vectors is localized in many complex networks 
these network models. 

We also consider here the case for which the adiabatic assumption is not valid. If the 
node degree is large enough, then x changes at a much slower rate than that of x. Therefore, 
this assumption will be valid if the degree of each node is large. On the other hand, if the 
degree of each node is small, then not only the adiabatic assumption but also the use of the 
central-limit theorem cannot be justified. For example, our theory will not work well for the 
"star" -network in which almost all nodes have a degree of 1. 

We note that theory for the BM model in the wealth-condensed phase is still elusive, 
as the standard central-limit theory is not applicable when (x 2 ) diverges. Thus, to discuss 
the static properties, we will need to generalize the central-limit theorem so that it can be 
adopted even if the variance diverges. 

Finally, we note that our developed technique is applicable to other dynamical systems 
subject to noise. Stochastic dynamics on complex network is attracting the interests of 
many reserchers recently, however, dynamics with multiplicative noise has not studied so 
far. We believe this method will shed new light on other problems regarding the dynamics 
on complex networks. 
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